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The Sun's global meridional circulation is evident as a slow poleward flow at 
its surface. This flow is observed to carry magnetic elements poleward - producing 
the Sun's polar magnetic fields as a key part of the 11-year sunspot cycle. Current 
theories for the sunspot cycle assume that this surface flow is part of a circulation 
which sinks inward at the poles and turns equatorward at depths below 100 Mm. 
Here we use the advection of the Sun's convection cells by the meridional flow to 
map the flow velocity in latitude and depth. Our measurements show the largest 
cells clearly moving equatorward at depths below 35 Mm - the base of the Sun's 
surface shear layer. This surprisingly shallow return flow indicates the need for 
substantial revisions to solar/stellar dynamo theory. 

Subject headings: Sun: dynamo, Sun: rotation, Sun: surface magnetism 



INTRODUCTION 



The Sun's meridional circulation has been a part of solar magnetic dynamo theory 
for half a century. A poleward meridional flow from the Sun's mid latitudes was invoked 
in the earliest models (long before the flow was actually measured) to transport magnetic 
elements from decaying sunspot regions to the poles where they would erode the opposite 
polarity magne tic field from th e old sunspot cycle and build up the polar fields of the new 
sunspot cycle (IBabcock 1 119611 ). This surface meridion al flow, along with its latitudinal 



structure and variation in t ime, is now well obseryed (ITopka et al. 



19931 : lHathaway et al. Ill996t lHathawav fc Rightmire 



magnetic flux transport is we ll established ( 



20K 



d. boiih 



DeVore fc Sheeley 



1982; Komm et al. 



and its role in the surface 



19871 Ivan Ballegooijen et al. 



19981 : ISchrijver fc Title Il200ll ) . Dynamo theories over the last decade and a half have assumed 
that the mass traveling poleward in the surface layers sinks inward at the poles and returns 
to the equator along the base of the Sun's convection zone at a depth of 200 Mm. In these 
theories this slow, dense, equatorward flow is responsible for the equatorward drift of sunspot 



-2 - 



activity (IDikpati fc Choudhuri Ill994t iDikpati fc Charbonneau Ill999t iNandy fc Choudhuri 
200ll ) . Dynamo models based on this deep meridional circulation have recentl y been used to 



predict the currently emergin g sunspot cycle, albeit with disparate results (IDikpati et al. 



20061 : IChoudhuri et al. 1120071 ) 



The meridional flow is difficult to measure. Its amplitude (~ 10 — 20 m s _1 ) is more 
than an order of magnitude weaker than that of the other major flows observed on the sur- 
face of the Sun. The axisymmetric longitudinal flow, differential rotation, has a dynamic 
range of ~ 200 m s _1 and the non-axisymmetric cellular convection flows have typical ve- 
locities of several hundred m s _1 . The meridional flow is observed to be poleward from 
the equator with peak flow speeds in the mid latitudes. The flow a mplitude measured from 
Doppler shift s of spectral lines formed at the surface is ~ 20 m s" 1 (IHathawav et al. Ill996l 



Ulrich Il2010l ) while th e amplitude found by measuring the motions of the sm all magnetic 
features is ~ 12 m s _1 (IKomm et al. Hl993l : lHathawav fc Rightmire ll2010L l201lh . The merid- 
ional flow can also be measured using the methods of local helioseismology which yield a 
peak velocity of ~ 20 m s" 1 that appea rs to be constant with depth down to ~ 26 Mm 



(IGiles et al. II1997I ; ISchou fc Bogartlll998l ). Recently, however, two new measurement meth- 
ods have indicated a decrease in ampli tude with depth. A method using global helioseismol- 
ogy (IMitra-Kraev fc Thompson 1 120071 ) found a meridional flow that decreased with depth 
and became equatorward at a depth of only 40 Mm - but with a large range of error. A 
method using the movement of the larger solar convection cells, supergranules, also found a 
meridional flow that decre ased with depth but with out precise depth information and without 
detection of a return flow (IHathaway et al. II2010I ). While some helioseismic studies i ndicate 
a poleward meridional flow at depths well below 26 M m. iDuvall fc Hanasoge I (120091) found 
that those methods are prone to systematic errors and Gough fc Hinc man ( 201oh conclude 



that the flow below 30 Mm remains unknown. Furthermore Beckers 



(120071 ) has suggested 



that projection effects may have also compromised some of the local helioseismology results 
and concludes that the meridional flow velocity may decrease with depth. 

Here we measure the meridional flow by tracking the motions of supergranules, but ex- 
tend the analysis to include larger cells with deeper roots. Numerical models of comp ressible 
convection with radiative transfer in the near surface layers (jStein &: Nordlund 1120001 ) clearly 
show that small cells dominate at the surface and larger structures are found at increasing 
depth. Supergranules cover the surface of the Sun and have a broad range of sizes that sam- 
ple a corresponding range of depths. We measure the motion of the pa ttern of supergranules 
by analyzing data acquired by the Michelson Doppler Imager (MDI) (IScherrer et al. Ill995[ ) 
on the ESA/NASA Solar and Heliospheric Observatory (SOHO) satellite in 1996 and 1997. 
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DATA PREPARATION 



The data consist of 1024x1024 pixel images of the line-of-sight velocity determined from 
the Doppler shift of a spectral line due to the trace element nickel in the solar atmosphere. 
The images are acquired at a 1 min cadence. We average them over 31 min with a Gaussian 
weighting function which filters out any velocity components that vary on time scales less 
than about 16 min. We then map these temporally filtered images onto a 1024x1024 grid 
in heliographic latitude from pole to pole and in longitude ±90° from the central meridian 
(Figure 1). This mapping accounts for the position angle of the Sun's rotation axis relative 
to the imaging CCD and the tilt angle of the Sun's rotation axis toward or away from the 
spacecraft. Both of these angles include mo dification for the most recent determinations 
of th e orientation of the Sun's rotation axis fjBeck & Giles I l2005t lHathaway fc Rightmire 
2010l ). We analyze data obtained during two 60+ day periods of continuous coverage - one 
in 1996 from May 24 to July 24 and another in 1997 from April 14 to June 18. 

We also generate and analyze simulated data to assist in our determination of the 
representative depths. We construct the simulated data from an evolving spectrum of vector 
spherical harmonics in such a manner as to reproduce the spa tial, spectral, and temporal 
behavior of the observed cellular flows (lHathaway et al. 1 120101 ) . The cells are advected in 
longitude by differential rotation and in latitude by meridional flow, both of which vary with 
cell size. 



3. DATA ANALYSIS AND RESULTS 



We determine the motions of the cellular patterns in longitude and latitude for strips of 
data by finding the displacement of the maximum in the cross-correlation with similar strips 
from images acquired at time lags of 2, 4, 8, 16, 24, and 32 hr. Each strip is 11 pixels or 
~ 2° high in latitude and 600 pixels or ~ 105° long in longitude. We repeat this procedure 
for 796 latitude positions between ±70° latitude and for each hour over the 60 days of each 
MDI dataset and 60 days of simulated da ta. This cross- correlation method was first used to 
determine the equatorial r otation rate by iDuvall I ( )1980l ) who concluded that larger cells live 



longer and rotated faster. iBeck & Schou I ( 120001 ) used a 2D Fourier transform method and 
found a rotation rate that increased with the wavelength of the features. 

We calculate the average differential rotation and meridional flow profiles and fit them 
with 4th order polynomials in sin# where 9 is the heliographic latitude (Figure 2). The 
rotational velocity increases with increasing time lags to a maximum at 24 hr but then 
decreases at 32 hr. The meridional flow velocity decreases with time lag and, at the 32 hr 



-4- 















J 








- 











Fig. 1. — Heliographic map details of the line-of-sight (Doppler) velocity from SOHO/MDI 
(top) and from the data simulation (bottom). Each map detail extends 90° in longitude 
from the central meridian on the left and about 35° in latitude from the equator (the thick 
horizontal line). The mottled pattern is the Doppler signal (blue for blue shifts and red 
for red shifts) due to the supergranule convection cells. The latitudinal movement of these 
supergranules yields a measurement of the Sun's meridional circulation. 

time lag, reverses sign. The cellular flows that live long enough to be positively correlated 
32 hrs later are moving equatorward. This is a clear detection of the meridional return flow. 
The individual data points have a standard error of ~ 1 m s _1 but the vast majority of 
points indicate an equatorward flow significantly bigger than this. The curves fit through 
the data points indicate an equatorward return flow of 1.8 m s _1 with a standard error of 
< 0.2 m b -1 . 

We determine the characteristic convection cell wavelengths for the different time lags 
using the wavelength dependence of the differential rotation and meridional flow profiles 
used in the simulation. The differential rotation measurements are largely reproduced with 
a relatively simple rotation velocity, u, relative to an inertial (sidereal) frame of reference 
given by 
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u(6,\) = [(1980 - 246 sin 2 - 365 sin 4 9) cos 9} 

[1.0 + 0.046 tanh(A/35)] m s" 1 (1) 

where the wavelength, A, is given in Mm. (Note that the prograde velocities plotted in Figure 
2 are relative to a frame of reference rotating at the Carrington rotation rate with uc{9) = 
1991 cos 9.) The meridional flow measurements are largely reproduced with a northward 
velocity, v, given by 



v(6,\) = [(65 sin 9 - 78 sin 3 9) cos 9} 

[tanh((35 - A)/20)] m s" 1 (2) 

where the reproduction of the return flow is particularly sensitive to the zero crossing occur- 
ring at ~ 35 Mm. 

Figure 3 shows the equatorial differential rotation velocity relative to the surface and the 
meridional flow velocity at 30° latitude along with the data points from the MDI (averaged for 
the two years) and simulation data analyses. The data points from the simulation virtually 
coincide with those from the MDI data except at the 32 hr time lag and for the northward 
velocity at the 16 hr time lag. No doubt, better fits could be obtained with more complicated 
flow profiles. It is apparent, however, that the differential rotation velocity must decrease 
for wavelengths greater than ~ 35 Mm and that the reversal in the meridional flow direction 
must be more abrupt. 

Figure 3 also shows that the points do not fall right on the curves for the input flow 
profiles. We attribute this to two different processes. The two latitudes chosen for this 
figure represent the latitudes at which each flow reaches its maximum. Since the convection 
cells span a finite range of latitudes the measured values should be less than these maxima. 
Howe ver, the differential rot ation signal is subject to an additional line-of-sight projection 
effect ( Hathaway et al.lbood ) which makes the Doppler features appear to rotate faster. This 



effect raises the measured values to increasingly higher values with increasing wavelength for 
the prograde velocity. 



4. CONCLUSIONS 



The relationship between the wavelength of a convection cell and the depth at which 
it is anchored or steered is well constrained by the stability of the surface shear layer and 
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observations of the rotation rate with depth from global helioiseismology. An increase in 
rotation rate with depth has long been suggested by observations and is attributed to the 
conservation of angular momentum for fluid elements moving inward and outward in the near 
surface layers (IFoukal fc Jokipii Ill975t iGilman fc Foukal Ill979l ; lHathaway Ill982l ). However, 
a rotation rate which increases inwar d faster than that give n by the conservation of angular 



momentum is dynamically unstable ( Ch andrasekhar 



rate increase from helioseismology (ISchou et al. 



1998 



1961). Measurements of this rotation 



Corbard fc Thompson 1120021 ) indicate 



that it follows this critical gradient to depths of 10-15 Mm. This gradient (the dashed line 
in Figure 3) and the helioseismicly determined rotation rates (open circles in Figure 3) are 
matched by the simulation input rotation profile if we assume that the convection cells are 
anchored at depths equal to their widths. Shallower cells would give unstable gradients. The 
mass density increases nearly quadratically with depth so it is reasonable to expect the cells 
to be advected by the flows near their deepest extent. 

This association between cell wavelength and depth indicates that the poleward merid- 
ional flow seen at the surface reverses at a depth of 35 Mm - the base of the surface shear layer 



where the rotation rate reaches its maximum. Although this s 



rallow return flow violates the 



assumptions of flux transport dynamos ( 


Dikpati & Choi) 


idhuri 


1994 


Dikpati & Charbonneau 


1999: Nandv & Choudhuri 2001: DikDati et al. 


20061: 


Choudhuri et al. 


20071). it was pre- 


dieted by numerical simulations of the effects of rotation on supergranules 


( Hathawav 


1982). 



it is in agreement with global helioseismology, and it helps to reconcile other observations. 

The surface has the slowest rotation and the fastest meridional flow. Small magnetic 
elements rotate fa ster than the surface and have poleward meridion a l flow which is slower 
than the surface ( iKomm et al. I Il993l ; lHathaway fc Rightmire Il2010l . l201ll ). Both velocity 
components for the small magnetic elements are matched at a depth of about 15 Mm. While 
supergranules do have a broa d range of sizes, their spectrum exhibits an excess of power 
at wavelengths of 30-35 Mm (lHathaway et al. I l2000h . Our results indicate that cells this 
size have depths roughly equal to the depth of the surface shear layer. This hardly seems 
coincidental but rather suggests an intimate connection between the characteristic size of 
supergranules and the depth of the surface shear layer. 

A meridional circulation confined to the surface shear layer would also explain why 
numerical simulations of the solar conve ction zone below th is surface shear layer have had 
difficulty producing the observed flows (IMiesch et al. Il2000l ). In particular, the meridional 
circulations in these simulations are highly structured in latitude and highly variable in time. 
The source of this structure and variability can be attributed to the small number (~ 100) 
of convection cells that populate the simulated volume. A meridional circulation driven by 
(~ 10000) supergranules - the convection cells that populate the surface shear layer - is far 
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more likely to be less structured and variable. 

This detection of a shallow equatorward return flow for the Sun's meridional circulation 
indicates the need for a reassessment of solar dynamo theory. The flux transport dynamo 
models, all of which assume and require a deep meridional flow, apparently cannot be correct. 
While other dynamo models exist, the majority of these place the dynamo action at the base 
of the convection zone with a possible secondary and less organized dynamo in the surface 
shear layer. Upon reassessment we may find that the surface shear layer plays a far more 
important role in the global dynamo. 

The author would like to thank NASA for its support of this research through a grant 
from the Heliophysics Causes and Consequences of the Minimum of Solar Cycle 23/24 Pro- 
gram to NASA Marshall Space Flight Center. He is also indebted to Ron Moore and Lisa 
Rightmire who read and commented on the manuscript. SOHO, is a project of international 
cooperation between ESA and NASA. 
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Fig. 2. — Flow profiles as functions of latitude from the 1996 SOHO/MDI data (top row), 
the 1997 SOHO/MDI data (middle row), and the simulated data (bottom row). The flow 
velocities measured at each latitude are shown with colored dots for each time lag as indicated 
in the figure. The solid lines with the same color coding represent the 4th order polynomial 
fits to each profile. The meridional flow decreases in amplitude with increasing time lag and 
reverses direction for 32 hr lags. The rotation rate increases as the time lag increases up to 
24 hr then drops at 32 hr. 
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Fig. 3. — The simulation meridional flow speed at 30° latitude (top) and equatorial differ- 
ential rotation relative to the surface (bottom) as functions of wavelength are shown by the 
solid lines. The observed values from the cross-correlation analysis with the MDI data and 
the simulated data (filled circles and crosses respectively) are shown at their characteristic 
wavelengths. Error bars centered on each symbol represent 2a errors. The dashed line shows 
the theoretical limit to the rate of increase in rotation rate at the equator and the open circles 



show the rotation velocity determined from global helioseismology by ISchou et al. I (1l998l ) 
both assuming that the wavelength equals the anchoring depth of the cells. 



